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SUMMARY 

The design of a large high-speed airship is primarily 
a structural problem, in which the most important 
stresses are those due, directly or indirectly, to aero- 
dynamic forces on the surface of the hull. The force 
on any small element of the surface is most conveniently 
divided into two components, respectively tangent and, 
normal to the surface. The tangent or skin-friction 
forces are so small- per unit area that they are structurally 
almost negligible compared with the normal forces; yet 
their total integrated resultant is responsible for almost 
the entire drag of the hull, whereas the normal components 
of pressure are so nearly balanced over a good hull that 
their net resultant is practically zero. The interreadion 
of these very substantial forces is, of course, through the 
medium of stresses in the hull, and in combination with 
iin and inertia forces they are essential not only from a 
structural standpoint but also in the consideration of 
stability and control. Tlie distribution of velocity and 
skin friction can also be indirectly determined from the 
normal force distribution. An accurate determination of 
the latter and its effects is therefore of the very first 
importance. 

The pressure on ellipsoidal shapes is 'presented first in 
Part I as a foundation for more generalized formulas. 
Although any ellipsoid is susceptible of accurate mathe- 
matical treatment, only the case of a prolate spheroid 
with circular cross section is investigated here because 
of its approximation to airship hulls. 

Part II deals with important adaptations of the ellip- 
soidal formulas, and other hydrodynamic relations to any 
airship hull, with particular reference to structural 
requirements. 

In Part III the theoretical results are applied to the 
practical computation of airship stability, and relations 
established which can be evaluated from simple wind- 
tunnel tests. 

In Part TV the same fundamentals are used in the 
determination of viscous forces, leading to an improved 
classification of airship drag, and a new outlook on drag 
generally. 

Examples of practical airship characteristics are em- 
ployed throughout. 



INTRODUCTION 

A determination of pressure distribution in a wind 
tunnel usually is tedious, expensive, and inaccurate. 
Even full-scale pressure readings have so far failed to 
give very consistent results, presumably due to local 
interference effects. Enough has been done experi- 
mentally, however, to show conclusively that the normal 
pressure is substantially unaffected by viscosity and 
skin friction except near the stern. The same con- 
clusion is also reached from considering the above 
mentioned small magnitude of the unit tangential 
force. In most of the work the viscosity may there- 
fore be ignored; and this is exactly the assumption on 
which the study of classical hydrodynamics has been 
based. This time-honored science, which for centuries 
was not much more than a mathematical toy, has thus 
found for itself at last a directly practical outlet in its 
increasing aeronautical use, particularly for lighter- 
than-air. 

The aerodynamic forces on airship hulls have lately 
been investigated by different methods based on 
hydrodynamical theory of flow. Professor Von Ear- 
man (reference 5) has applied the method of sources 
and sinks combined with assumed vortices at the stern 
to the investigation of pressures on model of airship 
Los Angeles. Dr. R. Jones (reference 3) determined 
the pressure distribution on a prolate spheroid based on 
classical hydrodynamical theories. Dr. M. Munk 
(reference 2) by applying the theory of momentums has 
derived the resultant moment acting on the airship; 
and by assuming 2-dimensional transverse flow, has 
approximated the distribution of transverse force 
acting on airship hulls. Some of these investigations 
have been presented in complicated manner and some, 
by disregarding certain factors, lead to erroneous 
conclusions. 

This paper, submitted at the request of the National 
Advisory Committee for Aeronautics, aims to over- 
come certain difficulties and to bring the subject gen- 
erally up to date. The purpose of the first two parts 
is to present in concise shape all the formulas required 
for computation of the hydrodynamic forces, so that 
they can be easily computed for either straight or 
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curvilinear flight. Improved approximations are also 
introduced having a high degree of accuracy through- 
out the entire range of practical proportions. The 
remaining two parts of the report are devoted respec- 
tively to stability and skin friction, as functions of the 
same hydrodynamic forces. 

SYMBOLS USED 

P= any point on the surface, 
P 0 = point at zero velocity, 

p = density of air (standard =0.00237 slugs), 

yx =■ coefficient of viscosity, 
V= velocity at point P (feet per second), 
V 0 = velocity of hull relative to air at oo (feet per second), 

2»= pressure at point P (pounds per square foot), 

2 = 2W = pressure at P 0 = | V Q 3 (pounds per square 
foot.) 
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■fineness ratio 
Fiqtoe L— Inertia coefficients of ellipsoids 

Coefficients of additional mass of ellipsoids: 
e 3 (w 2 -!) 3 * 



6 7i 2 V^- l-7ll0g(»+ V^-l) 

e 



where, 
e= eccentricity of ellipse 

n = fineness ratio = afb, 
a = major axis of ellipsoid, 
6 = minor axis of ellipsoid, 

B= l+kn, and is such quantity that^+^j^l, 

r -7S na+1 i T , 2 ^fr»+l)-6 
° K n'-l~ i ~ i 4An 3 -3(n 3 +l)' 

k i} k 3 , and k' may be taken from tables of Lamb and 
Munk, for A, B, and O-see Table I and Figure 1, 
0= angle the direction of wind makes with axis 
(pitched flight), 



^= angle the direction of wind makes with axis 

(yawed Alight), 
i^o = angle the direction of wind makes with axis 

(circular flight) — (at center of volume), 
a =» angle the tangent at P makes with the major 

axis, 

ao= angle the tangent at P 0 makes with the major 
axis, 

0= angle of arch between point P and a meridian 
perpendicular to the plane of the turning' 
circle, or the horizontal meridian (equator) 
in case of pitched flight, 
x,y = profile coordinates of the point P with center 
of ellipsoid taken as origin. (Fig. 2), x posi- 
tive backward, 

r = y = radius of circular cross section at point P, 
F 0 =P« (for curved flight), where «= angular ve- 
locity of ellipsoid (radians per second), 
P= radius of circle in which center of vol- 
ume moves (R=hull drag in Part IV), 



R 



/4x/S A- A is perpendicular 
to R and axis of ship 



A— 





Fiotoe 2.— Geometrical notation 

S = cross section area of ellipsoid at point P^irr 3 , 
S mK -t = maximum cross section area of ellipsoid = 
<LSycLc=rate of change of cross section =2 irr tan a, 
vol = volume of the hull, 
AF T = transverse force per foot of axis, 
AP £ = longitudinal force per foot of axis, 
AMi= longitudinal moment per foot of axis, 
F L = longitudinal force from bow to any station x, 
M=M£+ilf r =moment about station a; of all forces 

forward of that point, 
M 0 = total tanming moment about center of volume, 
M T = moment about station x of transverse forces 

forward of that point, 
M £ = same, of longitudinal forces, 
Q = transverse shear from bow back to any sta- 
tion x. 

See Figure 10 for additional symbols applying to 
Part HI. 
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PART I 

HYDRODYNAMIC FORCES ACTING ON A PROLATE SPHEROID 



For the fundamentals on which this study is based 
it is suggested that the references, particularly ref- 
erence 1, be consulted. It is sufficient to state here 
that the velocity distribution assumes frictionless, 
incompressible, irrotational flow (whether the flight 
path is curved or straight), as given by Laplace's 
equation for velocity potential $>. The pressure at 
any point P is then given by the expression: 

p d$ r? 
p^bt 2 

where v is the air velocity at P relative to undisturbed 
air. This is equivalent to an adaptation of Bernoulli's 
equation: 

whore U is the free air velocity of point P (=y 0 for 
a point at radius B from center of turn), or for any 
point on the hull in straight flight. 

NORMAL PRESSURE DISTRIBUTION 

For curvilinear flight, the most general case (see 
reference 3), the pressure at any point can be alge- 
braically expressed as: 

£=^cos ^o-jjsin ^ +Q^+sin^ 0 ) — |Acos^ 0 cosa 

+5sm^ 0 sin a sin <j>+[2S(C- V-CSn J** " Sm 

-(B^ 0 +^J cos 3 4, (1) 

(See Pt. in for the relation between f 0 and B.) 

In case the ellipsoid is moving straight at a certain 
angle of pitch 0, B = <» , fa = e, and the expression reduces 
to: 

^ = 1 - B 3 sin 3 0 cos 8 <f> — (A cos 0 cos a + B sin 0 sin a sin <£) 2 

In case 0=0, the expression for pressure at any point 
reduces to: 

2=1 -A 2 cos 2 a and%-*=A (2) 
Point P 0 , where velocity = 0, and pressure is a maxi- 



mum 



^2= -will be located where cot <*o = ^[ * an 0. 



At the point of Tmrnmum pressure where a=^ — ao 
1-A 3 cos 2 0-B 2 sin 2 0 



2 

= 1- 



A 2 cos 2 0 



1- 



Bhm 2 0 



sin* ao cos* ao 

At point of zero pressure, where V= V 0 
V 3 

y%* &m 3 (oq— a) = l 

when 0 = 0, the zero pressure occurs at the point where 

1 

cos a= ~^ 

The pressure at any point for pitched flight in the 
plane of symmetry may also be expressed as: 

1 



y^-sha 2 (ao— a) 



If we designate the side of ellipsoid turned to the 
flow as "windward," and the pressure at this side 
p w ; the side turned away from flow "leeward," and the 
pressure p t ; calling the meridian plane between the 
sides equatorial and expressing pressure on meridian 
plane p,; then in case of pitched flight: 



^=1— (4 cos a cos 6+B sin a sin 6) 3 



2 



• 1 — (A cos a cos 0— B sin a sin 0) 2 

^=1- J B 2 sin 2 0-A' cos 2 a cos 2 0 
2 



(3) 



(3a| 



(3b) 



If the values of p a , p t , and p, are plotted along the 
diameter of symmetry they do not he on a straight line, 
except for a very elongated shape. (See fig. 3.) 

At the equator the flow is at an angle ib' to the 
meridian such that, 

,, B tan 0 

tan = n 

v A cos a 



and the resultant velocity ratio, 
V* 

7 i 



frz^B 3 sin 2 6+ A 3 cos 2 a cos 2 0 
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.INTEGRATED FORCES 

Consider an element on the surface of an ellipsoid of 
revolution. The force due to aerodynamic pressure on 
a small arc per foot of axis mil be dF=pr d <f> 

This force can be subdivided into three compon- 
ents (fig. 2) : 

(a) Transverse, perpendicular to the axis AA. 

dF T =pr sin <£ d <£ 

(6) Longitudinal, parallel to the axis of ellipsoid. 

dF L —pr tan a d 4> 

(c) Perpendicular to plane of symmetry, parallel to 
the axis AA; this component is evidently balanced 
by the component on the other side of the hull and is 
here disregarded. 




(4) 



-SO -20 -IX) O +1-0 +2X3 +3.0 +4.0 

Pressure. Sx/sq. ft. 
Fiottbe 3.— Hydrodynamlc pre l im ) distribution over two cross sections— 
ZMC-S plotted against the diameter of symmetry; 9-10°,s-60m.p.h. 

The resultant transverse force per foot of axis on 
any section perpendicular to the axis, in case of circular 
flight is, 

j| - AF T ~2j^pr sin tfdtf = 2 [aB Jjf sin fa 

- ^ [(2AC- 2A+ sec 2 a) S- AOS^ cos 2 a cos fa 
In case of pitched flight, 

E = co, 

fa=0, 

and, 

AF r =2-2-^cos 2 ct sin2e=g — irrs]n2flsui2a (4a) 

The longitudinal component of pressure produces a 
longitudinal force per foot of axis, 



AF, 



= 2 pr tan ad<j>. 



In case of circular flight this evaluates to : 
AJf £ = 2 {2 cos 2 fa + 2 (g+ sin faJ+(£J 

— 2 (A cos 1^0 cos a) 2 — ^Bsin^o+^C^ -|~B sin fa sin a 



+ (20S- OS^- 2S) X ^JJtt tan a 



(5) 



In case of pitched flight the same expression re- 
duces to: 

AF L = q{2-2A 2 cos 2 6 cos 2 a 



— 2?(l + sin 2 a) sin 2 0} ut tan a 
In case of straight flight at 0=0, 



(6a) 



AS 



AF L =2q(l-A 3 cos 2 a) tit tan a = p (6) 

The effect of this force, although often disregarded in 
the present design of airships will undoubtedly be con- 
sidered with increase of speed. 

The longitudinal force per foot of axis if integrated 
over the length of ellipsoid will give the total longitu- 
dinal aerodynamic force to any section. If integrated 
over the whole length of the hull it will be equal to : 



Fx, => 2s (ki) (jjij^ sin fa, 



in case of circular flight. In case of pitched flight it 
will return to zero at the end. 

The total longitudinal force will reach its maximum 
compressive value a short distance behind the bow. 
At the maximum section it has a slight negative (ten- 
sion) value, which in case of straight flight of the ellip- 
soid (0=0) is, 

■ [l-^(l-2|^)J 

The expression of transverse force if integrated 
either graphically or mathematically over the length 
of ellipsoid from the bow to any station x produces the 
curve of aerodynamic shear. At the maximum sec- 
tion the shear reaches the value: 

Q~<L- 2 -S nax ^ ri {l-2 ¥Z ^} 

sin 2^o—22 h cos fa 

By integrating the shear curve over the total length 
the expression for aerodynamic moment around the 
stern due to transverse force can be obtained. 
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71? (Vol) 

M T - 2 (fa - fa) (vol) tfZT[ sin 2 to - 2qa ^ &i cos to 

The last two expressions for shear and transverse 
moment are derived for the case of circular flight; 
in case of pitched flight (to"=0)j the last terms of these 
equations will evidently be equal to zero. 

The longitudinal component of pressure, neglected 
in some of the recent investigations, produces a reverse 
moment, dM L =r sin $ dF L =pr 3 tan a sin 4> d <j>=r 
tan a dF T . 

Then the longitudinal moment per foot of axis, 
r tan a dF T = — rA F T tan a 

-x/2 

This moment will always have a sign opposite to the 
sign of transverse shear. If integrated over the length 
of the ellipsoid for either pitched or curved flight, it 
produces a total longitudinal moment, 

Mi=2(fe-fci) $-1) S" 1 2 

Combining this moment with the transverse mo- 
ment, the expression for total turning moment around 
the stern to which the ship is subjected, 

M=g.(fa-fa) (vol) sin 2 to _ 2 g_ fa (vol) ^ cos to- 

If the moment is taken around the center of volume, 
the last term disappears, making the total moment for 
either circular or pitched flight: 

M 0 =q(.fa-fa) (vol) sin 2 to- (7) 

This agrees with the expression derived by Doctor 
Munk in his' consideration of general flow around air- 



ship hulls. (Reference 2.) In his other work Doctor 
Munk has disregarded the effect of longitudinal com- 
ponents of pressure and, as will be shown, his expres- 
sion for transverse force is needlessly inexact. 

The fundamental difference between the circular and 
pitched flight in an ideal fluid is that in case of pitched 
flight there are no resultant transverse or longitudinal 
forces; while in case of circular flight there is a net 
transverse force component: 

F T =2g(Jc 1 ) (vol) ^L*° (8) 

and also a net longitudinal force component: 

2^=22(6,) (vol) (9) 

The two components have a resultant passing 
through the center of flight path curvature, thus 
satisfying the total energy conditions. 

It will be interesting to note that some of the aero- 
dynamic loads may be expressed as functions of 
"windward," "leeward," and "equatorial" pressures 
which were previously discussed, e. g., 

AF T =ir ^ 

AMl^h-^ (jpi—Jpw) tan a 

The above expressions are theoretically correct for 
circular as well as pitched flight. 
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APPROXIMATION OF HYDRODYNAMIC FORCES ON AIRSHIP HULLS 



A critical study of the mathematical formulas apply- 
ing to true ellipsoids throws considerable light on the 
effect of other curves. 

It may be observed that all the pressures and inte- 
grated forces are functions of the radius of cross section 
and the slopes at the particular point, and of general 
nondimensional coefficients of additional mass. Fur- 
thermore, there is a direct connection between these 
characteristics. Thus, for ellipsoids, it may be shown 
that 



J cos 1 a dx 

We now proceed to apply the expressions derived for 
an ellipsoidal shape to any airship hull. Several 
methods can be used for this purpose. 

Equivalent ellipsoids. — As the bows of many air- 
ships approach or are actually sections of ellipsoids 
(ZMC-2, R-101), these may be computed as ellipsoids 
and a modified curve can then be drawn for the stern 
of the ship. This method is recommended by Doctor 
Cox (reference 8), and will be discussed with the trans- 
verse stern force distribution. 

Another method for shapes departing more from the 
ellipsoid is to determine the additional mass coefficients 
for the complete hull and use them as constants in the 
ellipsoidal formulas, determining the actual airship 
forces from known characteristics of the airship curve 
(radius and slope). For the determination of addi- 
tional mass coefficients three simple methods of figuring 
equivalent ellipsoids may be used: 

1. Actual ft=2^ 

2. By length n-^/g^ 

3. By diameter n ^~~jjir 

The truth will usually be best expressed by the 
formula which gives the lowest result. Fine pointed 
shapes, for example, take the last formula, because the 
extreme point of the tail simply adds to the geometrical 
length without appreciably affecting the aerodynamic 
characteristics. The effective value of n is also de- 
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creased by the dissymmetry between the bow and 
stern profiles, but not as far as would be obtained by 
assuming the bow profile at both ends. 

A still more general, but still approximate method, 
is to fit at any point on the hull profile an ellipsoidal 
element, which would reproduce all the hydrodynamic 
characteristics at that point. This ellipsoidal element 
will, as far as possible, have the same geometric char- 
acteristics as the point on the curve, but tentatively, 
may be a part of an ellipsoid whose over-all dimensions 
do not correspond at all to the hull profile, i. e., its 
additional mass coefficients may be entirely different 
from the hull shape. The application of these equiv- 
alent ellipsoids may be checked by applying the general 
hydrodynamic proposition for frictionless, nonviscous 
fluid that integration of axial and transverse pressure 
components over the surface of any stiearnline shape 
must equal zero for uniform straight motion. This 
method will be more particularly applied to the deter- 
mination of transverse force to be described later. 

The pressure distribution plotted according to this 
method is shown on Figure 4, curve 1. This curve was 
plotted assuming the equivalent ellipsoid fitted to any 
point on the hull to have the following characteristics 
relative to the original hull curve : 

1. Concentric about longitudinal axis, 

2. Maximum diameter equal, (Z?=26,), 

3. Vertical ordinates at the point equal, (r ="}/,), 

4. Product of first and second derivatives of the 
curves at the point equal, 

which lead to the following relation between the fine- 
ness ratio of an equivalent ellipsoid and the charac- 
teristics of the point on the hull curve: 




PRESSURE DISTRIBUTION BY SOURCE INTENSITY 

For still greater accuracy we must go back to more 
fundamental relations. The method of deterrnining a 
stream-line form from a predetermined source and hi'tiIt 
line is not new and was used in Naval Architecture and 
Airship Design. (References 4 and 6.) It was expected 
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that the forma of low resistance could be determined 
by this method and that these forms would have an 
advantage that, knowing the hull profile and the in- 
tensity of sources and sinks, the pressure distribution 
could be easily found. In most cases, however, a 
predetermined distribution of sources and sinks does 
not give a profile curve with any simple mathematical 
equation that can be conveniently manipulated. 
Therefore, in practice it seems preferable to predeter- 
mine the profile curve and find the source distribution 
to correspond, though recognizing the greater initial 
labor involved. 

For determination of the sink and source line shown 
on Figure 4, Taylor's method (reference 6) was used, 



characteristics. The same procedure was then applied 
to the actual airship hull. 

These equations after being solved for/(2), repro- 
duced the sink and source line corresponding to the 
hull shape. Knowing /(2) values, the horizontal 
velocity at the point was determined by the equation, 

then the final velocity at the point, 

T7* VO+V* 



cos a 



This value of the velocity was used in the pressure 
equation, 




70 80 
Feet 



too no 



130 140 150 



Fiodee i. — Comparison of pressure carve derived from sink and source line (S & 8) with curve obtained by fitting ellipsoids. 9— 0°,p— Mm.p.h, 

ZMO-S shape, EH carve with E stem 



except that the problem was the reverse one of plot- 
ting the curve of source and sink intensity for a 
known body. This was done by means of the equation 
(for any one transverse section) : 

Where/(2), the source intensity, varies as a function 
of x (here called 2 to distinguish it from the abscissa of 
the section under investigation), and is integrated over 
the length of the hull in each case. 

Thus a complete series of equations was derived, all 
these equations being interconnected by the known 
equation of the hull curve. The method was first 
checked by reproducing a full ellipsoid having known 



I" 1 



for determination of the pressure at a given point. 

Von Karman (reference 5) made an independent 
investigation of pressure distribution using a sink and 
source line to determine the flow in a longitudinal 
plane, and a system of double sources to determine 
the flow in a transverse plane. By writing general 
expressions for the flow, he determined the values of 
the velocity components in both planes, due to each 
system of flow separately. Knowing the values of 
these velocities, the pressure was determined by 
Bernoulli's equation. 
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The pressure distribution curves based on ellipsoidal 
formulas for straight flight are plotted in Figure 4 for 
the metal-clad ZMC-S hull: (1) by the method of 
fitting ellipsoids to points on the hull; (2) by applying 
the mass coefficients of the complete hull to equation 
(2); for comparison (3), the actual hydrodynamic 
pressure distribution by the method of sources and 
sinks is also shown in the same figure. 

TRANSVERSE FORCE EST PITCHED FLIGHT 

The main difficulty in applying ellipsoidal elements 
lies in the approximation of additional mass coefficients 
to fit not only the point where the forces are investi- 
gated, but also to fit the complete hull. 

In this connection an interesting observation was 
made that the product of inertia coefficients AxB for 



AJV=2^cos 2 a sin 2d = <[ir r sin 20 sin 2a (10) 

The above simple, single formula is all that is neces- 
sary for plotting this most important aerodynamic 
force with a high degree of accuracy as to general 
distribution and magnitude. (Fig. 5.) After this force 
is plotted it will be found, due to the above approxima- 
tion, that the positive force over the bow and negative 
force over the stern do not balance exactly. From a 
quantitative standpoint this discrepancy is unim- 
portant, but if further analysis is required on the 
basis of an ideal fluid, it may be desirable to have an 
exact balance. This can be easily done, at the same 
time making the force curve at the bow a still closer 
approximation, by simply moving the whole curve 
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Fiqube 5.— Hydrodynamic loads, shears, and moments— ZMC-3 hull. EH curve with E stern, 9=10°, p=W m. p. h. Total longitudinal force computod for 9»0° 



all usable fineness ratios is practically equal to two, 
as is shown in Table I. 

Returning to the ellipsoidal formulas, it will be 
noticed that the mass coefficients are represented 

exclusively by the factor in the expression (4a) 

for transverse force per foot of axis. Thus for ellip- 
soidal type of flow the point characteristics (i. e., 
ordinate and slope) of any sectional element determine 
the transverse force almost independently of charac- 
teristics elsewhere. We may assume that the same 
holds true for a series of gradually varying elements 
and hence for the hull curve itself. 

AJ3 

It would thus appear that the factor -g- can be 

neglected, for most practical purposes, over the 
entire range of usable fineness ratios. The general 
formula for transverse force in pitched flight may 
then be written, 



bodily up or down enough to make the positive and 
negative areas balance. 

The close approximation of the lateral force derived 
by the above method to the force measured by wind 
tunnel tests can be seen from Figures 6, 7, and 8, 
showing the distribution of transverse force on modern 
airship hulls. 

Von Karman (reference 5), extending his investi- 
gation of pressure distribution to determination of 
transverse forces has arrived at the following expres- 
sion, 

AF T = pvr (UxWj + u r w T ) . 
Where u and w are the velocities due to the longi- 
tudinal and transverse motion, respectively, the sub- 
letters x and r refer to the respective longitudinal 
and radial components at a point in the plane of sym- 
metry. 

Applying this expression to ellipsoidal analysis, wo 
get, 



APPLICATION OP PRACTICAL HYDRODYNAMICS TO AIRSHIP DESIGN 



131 



AB 

UxWx+UrWr^—^- Vg sin 2 a sin 26 



or, 



approximately, — -^- sin 2 a sin 2 6 (the same result). 

AERODYNAMIC SHEAR IN PITCHED FLIGHT 

After the curve of transverse force is plotted, it can 
be graphically integrated by means of a planimeter 
and the curve derived will represent the aerodynamic 
shear Q to any point, in pitched flight. The maximum 
shear, at a point where AF T = o, may be approximated 
by a similar method to what was used for transverse 
force, thus: 



where, 



Qmni = a-yjB-A S^t sin 20 



is the mean effective value using apparent mass coeffi- 
cients for the entire hull. 

LONGITUDINAL MOMENT IN PITCHED FLIGHT (PER 
FOOT OF AXIS) 

As was mentioned for ellipsoids, the longitudinal 
moment per foot of axis, 

AM L = — A F T r tan a 

Knowing therefore the value of AF r and the geo- 
metrical characteristics of the point, the expression 



the moment reaching a maximum at this point. (Fig. 
5.) The net value of moment curve ordinate at the 
stem represents the total hydrodynamic turning 
moment on the hull, and must satisfy the equation (7). 

M=g (fe-Jfe,) (vol) sin 2 6 




ga 0 



Fioube 6.— Transverse force on model BS-l (e-9°). L Transverse force from 
wind tunnel tests by Bureau of Standards (Technical Data Files, MoCook Field, 

DS27/BS-1). n=AJ>=sXg^|xsln 2 e cos> a (corrected for hydrodynamlo 
balance.) m=AJr=5X^ (*»— *i)Xsin 2 B Monk's formula (ij-iO-0.78 

TRANSVERSE STERN FORCE 

The plotting of hydrodynamic forces, shears, and 
moments does not correspond to actual conditions 
observed on the airship hull. It will be noticed that 
in an ideal fluid no resultant force occurs, and the air- 
ship is subjected only to the action of a couple, whose 
magnitude equals the hydrodynamic moment. Ac- 
tually, the area under the distributed transverse force 




Fioube 7.— Transverse force on model B-33 (9-10°). L Transverse force from wind tunnel tests (A. B. O.; B & M No. 801); 
n-Air=8X j^Xstn 2 8 cos» a; m-APi—a jj| (h-ti) sin 2 S Munk's formula (ii-ii=.924) 



for longitudinal moment can be easily plotted. The 
longitudinal moment should have a negative sign over 
the entire length of hull. 

NET AERODYNAMIC MOMENT IN PITCHED FLIGHT 

The net moment to any point will be the algebraic 
sum or difference of the areas under the shear and 
longitudinal moment curves. This assumes tempo- 
rarily that the hull is held as a cantilever at the stern, 



curve at the stern is considerably smaller than at the 
bow, producing a resultant force at the bow in the 
direction of inclination of the ship's axis. (Figs. 6, 7, 
and 8.) Superimposed on the hydrodynamic flow at 
the stern there is evidently an airfoil type of flow in 
a direction opposite to the ship's inclination, producing 
vortices as is the case with an airplane wing. 

Doctor Von Karman (reference 5) has computed the 
magnitude of these vortices and has arrived at results 
closely approaching the actual load conditions. 
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The following method, somewhat approaching that 
recommended by Doctor Cox (reference 8), can be 
more conveniently used. In this case it is desirable 
to obtain from wind tunnel tests at high Eeynolds 
Number the resultant force and moment on the bare 
airship hull under conditions similar to those analyzed. 
When these values are not available, they can be 
approximated from results of tests of airship models 
of shapes similar to the one analyzed. The curve of 
hydrodynamic transverse force can be then modified 
at the stern in such manner that its area would be 
reduced by the amount equal to the force determined 
in the wind tunnel. The shape of the curve should be 
modified to approach the curves of transverse force as 
determined in wind tunnels on hulls of similar general 
characteristics. This distribution of transverse force 
will then approach very closely to the actual condi- 



conditions. Therefore, the investigation of pitched 
flight alone will usually serve as a guide for the deter- 
mination of the longitudinal strength of airships, al- 
though circular flight should also be at least tentatively 
investigated, if there are any doubts as to the 
strength in this plane. The relation between B and 
&>in actual flight and their effect on stability is treated 
in Part m. 

CONCLUSION 
The investigation of the possibility of applying the 
hydrodynamic ellipsoidal formulas to airship hulls 
leads to the conclusion that these formulas can be 
divided into two groups: (1) Formulas which can be 
approximated to permit the determination of flow 
characteristics from geometrical properties of the 
individual point at which the flow is desired, and (2) 
formulas which include the general shape as a whole. 




Piouee 8.— Transverse force on model B-101 (9-10°). Scale of model^one-quarter; 
length of model=33.01 Inches 

tions, and if the forces at the control surfaces and the 
inertia or other balancing forces are also considered, 
the complete load diagram can be reproduced. 

CIRCULAR FLIGHT 

For the determination of forces in circular flight, the 
exact ellipsoidal expressions should be used, making 
approximation for the additional mass coefficients as 
outlined above. The distribution of transverse force 
over the bow for circular flight of the ZMC-S is shown 
in Figure 9. This distribution is plotted for the same 
speed, with the angle of yaw at c. g. the same as in the 
pitched flight investigation. It can easily be observed 
that the transverse force is smaller in the case of cir- 
cular flight than in the corresponding case of pitched 
flight. Considering also that the transverse force in 
circular flight is balanced by inertia force largely dis- 
tributed over the hull, while the force in pitched flight 
may be largely balanced by excess of weight concen- 
trated at the points of load application, it can clearly 
be seen that pitched flight will produce higher shears 
and moments than circular flight for otherwise similar 



Figure 9.— Circular flight. Hydrodynamlo loads, shears, and moments— Bow ZMO-3 
(*=10°, s-M m.p. h., JJ-600 feet) 

The formulas which include only the geometrical 
point characteristics, such as the formula for trans- 
verse force in pitched flight, leading to the determina- 
tion of shear and bending moment, can be applied to 
any airship hull, and for the most part are probably 
even more accurate than the results of actual tests. 

The formulas for pressure distribution and for longi- 
tudinal force in the case of pitched flight, and for all 
the aerodynamic forces in case of circular flight, are 
dependent on characteristics of the point and also on 
the hull shape as a whole. In this case the designer 
should be careful in choosing the coefficients of addi- 
tional mass corresponding to the designed hull. If 
these are properly applied results very closely ap- 
proaching those of true hydrodynamic flow can be 
obtained. 

In regard to the stern force, it should be noted that 
the hull force in the immediate neighborhood of the 
fins, particularly between and behind them, is further 
modified by the presence of the Atir themselves. It is 
proposed to deal with this more particularly in a sub- 
sequent report on fin design. 
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AIRSHIP STABILITY 



The consideration of airship stability in yaw is 
based on the principle that any tendency of the air- 
ship to yaw from its path may be considered as the 
tendency to swing into a circular path deviating from 
the original direction. The study of stability in yaw 
therefore simplifies itself to a study of the airship in 
circular flight and a consideration of the transverse 
forces and moments acting on it. 

In this connection it should be first noted that for 
motion of the airship deviating from the straight path 
the angle of attack is no longer fixed but will vary 
throughout the length of the airship, reaching the 
highest values at the stern. The wind-tunnel tests 
give only the moments at a constant angle of attack 
throughout the hull; but as the principal stabilizing 
forces on the airship act at the stern, the wind tunne 
results may be utilized providing proper allowance is 
made for the above-mentioned variation of yaw angle 
in actual flight. This may be done by the use of hy- 
drodynamic principles. 

Consider an airship of mass m=pX vol, traveling 
on a circular path, with angle of yaw = ^ 0 (relative to 
undisturbed air); velocity tangential to path=y 0 and 
radius =22 (all taken at the center of volume). The 
center of gravity is assumed to be coincident with 
the center of volume, fa and fa are additional mass 
coefficients of the airship, A and B the total virtual 
mass coefficients. (Fig. 10.) 

The airship is assumed to be subjected to the fol- 
lowing forces: 

1. Transverse component of centrifugal force due 



to the mass of the ship itself- 



' B 



cos f Q distributed 



around the hull' 



through the hull, with resultant acting at center of 
volume. 

2. Transverse component of virtual centrifugal 
force due to the additional mass of air set in motion 



-g— fa cos which appears as aero- 
dynamic pressure with resultant applied at center of 
volume. This is equivalent to the ellipsoidal equa- 
tion (8), but may also be derived from a more general 
consideration of momentum. (Reference 2.) 

3. Distributed aerodynamic force which does not 
produce any resultant force but produces an unbal- 



ancing moment, M 0 = (fa— fa) sin 2 \J/ 0 ; from 
equation (7). 

4. Stern force acting at a distance I from the center 
of gravity, whose magnitude to balance the above 
ffiV 2 

moment should be -^f-ifa—fa) sin 2 

\(2) 



0) 



(3) W 



(3) 



(')= -ff-cos.fo 
mV* 

% " angle of yaw at e.g. 
R " radius of turn of e.g. 
m - p (vol.) = airship mass 
Y " yaw force measured in 
funnel 

N " yaw moment measured 
in funnel 
M 0 " moment due to distribut- 
ed force (3) 
C " distance to point of zero 
yaw from e.g. 

ofni 



*» distance to point of stern 

force application from eg. 
■ true airspeed tangential to 
path at e.g. m ft./sec. 
/c,,ic 2 - additional apparent mass 

factors (IJ,(2), (3), (4) = trans- 
verse forces. 

Figure 10— Transverse forces on airship hull in circular flight 

Eelation between turning radius and angle of yaw. — 
As long as the airship is in a condition of steady turn,, 
the forces should balance or, from the above, 

^a+*i)cos* 0 =^(k-*i)sin2*o. (H> 



thus, 



~Ttr~ — i — 8111 +° 



A 

r 



B-A . , 

■ — ^ — sm \{/ 0 
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from which, 



B= 



I I 



(k 2 - kt) sin f 0 (k 3 - h) ifo 



57.3- 



57.3^2 
'MB- A) 



(12) 



(^o hi degrees). 

It may he noticed that the angle of yaw of any one 
point on the axis is a linear function of distance along 
the axis. Thus as ^ 0 is the angle of yaw at the c. g., 
let c be the distance forward of c. g. to point where the 
local angle of yaw is zero. Then, 



c=B sin ^o* 



Al 
B-A 



ISxpfiB the local angle of yaw at the point of applica- 
tion of the stern force (I feet behind the c. g.), 

. . c+l. . B. . 
tan tan y 0 =-^tan \f/ 0 



or, 



(13) 



from which we get as an alternative value of the turn- 
ing radius, 

The above considerations of an airship in a perfect 
fluid in a condition of neutral stability are only approxi- 
mated in practice, due to the distributed nature of the 
tail force, as already noted. In the case of a bare hull 
at a fixed angle of pitch or yaw, the distributed hydro- 
dynamic force at the stern is reduced while the dis- 
tributed force at the bow remains the same. Usually 
this difference is termed "dynamic lift on bare hull" of 
the airship. Its real center of action is, of course, for- 
ward of the c. g. For the purpose of the above analysis, 
however, it is considered not in the sense of a resultant 
force at the bow of the airship, but as a force at the 
stern acting in a direction to oppose the theoretical 
turning moment M 0 which would otherwise apply in 
the analysis of a perfect fluid. This force is the Yj, 
force which is measurable in wind-tunnel tests on a 
bare hull for a certain angle of yaw f or pitch 0. Its 
moment arm about the c. g. as measured in the wind 
tunnel is, 

M S-N t 



"When the fin surfaces are introduced they produce a 
force Y f and moment Y f l f . For the purpose of analy- 
sis we here assume Z ft =Z / =Z combining Y h and Y f into 
a single force Y, which is not far from actual results. 

An experimental check is found in a comparison of 
the radius of turn computed from formula (12) with 
that obtained in full scale tests. In the case of the 
Los Angeles the computed radii are 4 per cent off from 
the experimental, but are closer to the experimental 



values than the results derived by the hydrodynami- 
cally inexact formula used in N. A. 0. A. Technical 
Report No. 333. The full scale tests on the 0-7 (N. A. 
C. A. Technical Report No. 208) are within 2.3 per 
cent of the computed figure. A somewhat larger dis- 
crepancy in the case of the ZMC-2 will be considered 
later in connection with the car effect. 

Stability criteria. — Comparing actual flight with the 
theoretical, where the forces were in balance, it oan be 
seen that if the actual force at the stern is larger than 
the theoretical stern force opposing the turning move- 
ment, the airship will tend to return from the circular 
to its original straight path. 

Whether this condition is fulfilled may now be deter- 
mined from a simple static wind-tunnel test, preferably 
at high Reynolds Number. The only further assump- 
tion involved is that the stern force for a given airship 
speed is determined by the angle of yaw at the fins, 
regardless of what it is elsewhere. As will be shown 
later, this seems to be almost exact for the fin force 
proper and usually a fair approximation for the balance 
of the stern force. Thus, for stability computation, 
the wind-tunnel angle ^ must now be taken as equiva- 
lent to ipf, not to yj/ B , because from (13) 



, A . 



A t 



(13a) 



Referring back to the condition for neutral stability, 
equation (11) for small angles can now be written: 



B 



(B-A) 



57.3 



of which the first term is the actual centrifugal force, 
including that of the additional air mass, and the 
second term is the balancing stern force necessary to 
hold a constant radius of turn. If Y exceeds this force 
we have positive stability. Generally speaking then, 
we may define the degree or criterion of stability as 
the ratio of the actual stem force to the balanoing 
stern force. Then any value greater than unity is 
stable. Thus we get stability criterion No. 1 : 



S. O, 



28MYB 



(14) 



,y q(voi)A(B-A)t 

The same method of reasoning can be also applied 
to moments. The moment measured in the tunnel 
for a certain angle of yaw will be evidently: 

N=M 0 f -YA-Y f l f 

The moment due to the stern force alone is evidently, 

-Mo'-N-YA+YJftsYl (15) 

This moment should be larger than the theoretioal 
moment to have the ship stable, or 

(M 0 '-N)>mVo 3 (B-A) jj 
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Taking the ratio as before, we get stability criterion 
No. 2: 

A third stability criterion may be based on an exper- 
imental determination of the resisting yawing moment 
by the use of damping apparatus or a whirling arm- 
It seems obvious, however, that additional mechanical 
complications of this kind are not justified unless they 
serve to reproduce more closely the type of motion 
found in actual flight. For example, the oscillating 
type of damping apparatus would seem to be par- 
ticularly futile. Though any method by which the 
stern force or moment can be plotted against the local 
angle of attack at the stern gives a means of approx- 
imating the required ratio, the most important test 
condition is high Reynolds Number. 

The evaluation of a damping test is best expressed 
as a moment N D around the c. g. (computed from the 
test moment around the actual center used) for a 
given transverse component of stern speed v f at the 

position I behind the c. g. Then ^ corresponds to 

(in radians), and we get as stability criterion No. 3: 

^^lO- py.voi^-A) J ^ 

Tail arm length. — Proper use of the formula for 
$. O.y depends on getting the correct value for I. 
The fin arm l f may be taken as the distance from the 
c. g. of the ship back to the mean c. p. of the fins, the 
latter point for any one fin being found as follows: 
Find the centroid of the entire fin, including control 
surface; draw a line through the centroid, in the plane 
of the fin, parallel to the hull meridian at the front 
of the fin; from where this line cuts the leading edge, 
measure back along the same line one-fourth of its 
length from leading edge to trailing edge. For most 
models I can be taken as approximately equal to l f 
but is more accurately given by the expression: 

l= M^N^(j_ yol {B _ A)xl/ _ N ^ (18) 

In case this gives a value differing much from l f , the 
bare hull force and tail arm, Y h and l h may be similarly 
determined from a test without the fins, and the figures 
checked by equation (15). In special cases portions 
of the hull force may be considered quite separately 
from the fin force, in a similar manner to the car force, 
now to be considered. 

Effect of car. — It is clear that an additional fin force 
at the center of gravity, such as due to a car, will serve 
to carry part of the centrifugal force. For a given 
angle of yaw it will cause the ship to turn on a smaller 
radius, which in turn will produce a larger local angle 
at the fins, thus improving the stability. Another, 
less precise way, of explaining it is to say that there is 
a relative damping moment between the car and the 
149900—33 10 



fins which opposes a tendency to turn. The important 
thing to note here is that the above two statements 
are not of two different effects, but refer to one and 
the same> effect. 

Though the more general case of a plurality of side 
forces at different points on the hull may be solved by a 
similar method, there seems to be no practical import- 
ance in trying to express the result in a single algebraic 
formula. Hence, we return to the case of a side force 
at the c. g. In this case, let Y c be the difference in 
transversive force with car on and off, measured at an 
angle of ^ in the tunnel. Then by the principles 
already laid down, the actual car force in free flight 

will be for an angle of ^ at the fins. The balance 

of forces for equilibrium is: 



where, 
and, 



(19) 



M 0 =mV 0 2 (B-A) 



57.(3 



From (19) we get, 



l m 28. QY c -4 B 
&>"*2vo1j4.^ A 

which can be substituted in the expression for 22. For 
a wind tunnel moment N, we have the second stability 
requirement as before: 



M-N> mV<? (JB—A) 



57.3 



(20) 



which, expressed as a ratio in terms of ^ is: 

Thus, aside from interference effects, the stability is 
slightly improved by a fin force at the c. g. Even 
further forward, there is not likely to be much negative 
effect, due to approaching the point of zero yaw. Put 
more broadly, it is on the safe side to consider any 
distributed load as concentrated at the point of appli- 
cation of its resultant. 

A word of warning here: measure N with the car (or 
cars) on, to be sure of catching any blanketing of the 
fins so caused. If Y is the stern force, S. O. K ' is 
derived from equation (20) by putting M— iV= Yl and 
substituting for ^ 0 in terms of i>. 

Pressure and velocity distribution. — Having deter- 
mined the free flight connection between f and B 
we may substitute either way in equation (1), Part I. 
A practical example is plotted in Figure 9. Expressed 
in terms of the stern angle if/, where x=lby means of 
equations (12a) and (13a), we have a convenient means 
of checking the validity of the assumptions on which 
the stability criteria are based. 
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The particular question here is whether there is 
enough difference, in the hitherto unconsidered, flow 
pattern, as between turning and simple yaw to make 
a serious difference in the resulting stern force for a 
given value of f. Taking the extreme condition of a 
fin element subjected to the full theoretical flow at the 
hull surface, its actual reaction, compared to the free 
air value, is increased in the ratio, 



where, ip' is the angle between the streamline and the 
meridian at the point in question. But for small 
angles (in radians), 

*'=sin*'=^ 

or, 

W e V m V c 

where V m is the meridian component and V 0 the cir- 
cumferential component of V. 

Equation (1), by its derivation, consists of four main 
terms of which the first two are the component terms 



of 



the third is 



and the fourth is 



V* 

w 



For the 



position of maximum fin effectiveness 4> = 0 ; also putting 
sin to—fa (radians), cos ^ 0 =1 and substituting for B, 
\f/ 0 and x as above we get, 



F=A(A+C-4jp) 



cos a 



Referring now to Table I or to Figure 1, it may be 
seen that the factor A (a+ @~^£f^ takes a value of 

2.25 for 7i=l, 2.0 for 7i= oo, and 1.94 for 7i=5, or is 
practically equal to 2 for the entire range of usable 
fineness ratios. 

Applying similar approximations to equation (3b) 
for simple yaw (0=^), we get as the factor of increase, 

F=AB cos a»2 cos a^F' 

Therefore, within the range of hydrodynamic condi- 
tions, the fin force at small angles is close to a straight 
line function of f for either curved flight or simple - 
yaw, as originally assumed. The boundary layer and 
other modifications of flow, probably of small effect 
in relation to the present subject matter, will be taken 
up in a later report. 

limitation of assumptions. — It should be noted 
finally that the assumed straight fine variation of the 
hull and fin forces with yaw angle is modified for the 
hydrodynamic forces when ^ becomes so large that the 
relation, 



sin 2^=2 sin ^=2^ 

no longer holds within the desired order of accuracy. 
And the assumption for the fin force breaks down 
entirely at the burble point of the fin airfoils. This is 
particularly noticeable for fins of high aspect ratio, 
such as used on the ZMC-2, this ship showing a very 
definite tendency to "spin" in an extreme turn. 

Conclusions. — For most purposes, equations (14) 
and (16) are recommended, giving greatest credence 
to the one giving the lowest value of S. O. and designing 
so as to make it at least 1.0 and preferably 1.15. 
Serious discrepancies between the two criteria should, 
however, be investigated for effects of distributed force, 
interference, and burbling. The angle is preferably 
taken as half of the algebraic difference between equal 
positive and negative angles, not more than ±10°, 
T and 2V being similarly derived. For larger angles, 

substitute ^ for ^ and ^ for ^ plotting the values 

of S. O. against \p. 

Stability in pitch may be analyzed in exactly the 
same way, the principal difference being the negligible 
effect of the cars, and the static righting moment due to 
the vertical distance between the center of volume and 
the c. g. For the latter reason there is seldom any 
^difficulty about the stability in pitch if that in yaw is 
satisfactory. 

The ZMO-S, designed by these methods, indicated a 
slightly positive stability at small angles, a result well 
substantiated in practice. Examples of unstable air- 
ships were the Shenandoah, and the Army AO. There 
is no question about the latter by results from either 
wind tunnel or actual flight. The Shenandoah, how- 
ever, although admittedly unstable in flight, did not 
seem to have as bad a reputation as its low criterion 
would indicate. Its very long and generally poor shape 
gives a distribution of transverse hull force sufficient 
to be a substantial factor. But even taking this into 
account the criterion is still quite low. Hence, it seems 
probable that the results of experience with the airship 
were largely colored by the large radius of turn, and 
by the time element contributed by the large moment 
of inertia involved in any unstable deflection from the 
straight path. In this case the principal danger of an 
unstable condition remains in its effect of magnifying 
the forces due to a sudden gust. Thus a sound and 
careful analysis of stability is of the utmost importance, 
particularly in the case of large airships where the 
result from a control standpoint is likely to be masked 
by elements involved in the size as such. 

In any case the establishing of stability is but a 
preliminary to the determination of actual forces for 
specific conditions of flight. 
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FBICTIONAL FOECES 



Here we determine the approximate distribution and 
magnitude of the frictional drag on a large airship hull. 
Analysis will proceed on the following assumptions: 

1. Turbulent type of boundary layer, of negligible 
thickness compared to the hull dimensions. 

2. Flow otherwise frictionless and irrotational. 

3. Distribution of unit tangential force proportional 
to the square of the local air velocity relative to the 
surface. 

4. Magnitude of the total integrated force deter- 
mined with reference to the mean effective Reynolds 
Number of the hull as a whole. 

5. Axial motion of an ellipsoidal hull. 

Foroe and energy relations. — Let dBo be the hori- 
zontal component of skin friction resistance on a cir- 
cular element of surface of axial dimension ds. In the 
preliminary analysis we assume that this force varies 
exactly as the area, and as the square of the air velocity. 
Therefore, we may write, dBo^GpryV 2 da;, where O is 
a coefficient, and p is the density in slugs. From the 

ellipse equation, y=^-Ja?-2?, by hydrodynamics, 
Eq. (2) 

W-A 2 Vo* cos 2 a 

_ A'Vo a (fl a -3?) 

Then, for the entire ellipsoid, 

2G P irA 2 V 0 i f a (a 2 -a?) 3 * ^ 
11 J V-s*(l-I) 

where V 0 is the speed of the ship (feet per second) and 
A is the hydrodynamic relation, -^p^l +fa. By in- 
tegration between the prescribed limits, 



(1) 



With the proper value of O, to be determined, the 
above equation gives the integrated axial force con- 



tributed by the skin friction. That it is not the total 
drag due to the skin friction may now be made appar- 
ent by energy considerations. 

Let cLE be the energy per unit time absorbed by the 
same circular element of surface, along the meridian 
increment ds. Then if dF is the force along ds (inte- 
grated around the circle) ; 

d_B= VdF<= Cp-rj/F 3 di " 



cos a 



(2) 



but, 



V S =A 3 V<? cos 3 <*• 



AWf^-i?) cos a 




FlffUBE 11 

Substituting and integrating as before, 

E-^OpAWW^ 

But if the total drag is B, E= VoB, 
and, 

R^OpAWo'a'^jy (3) 

or, 

In other words, the skin friction drag B 0 must be 
increased by a "pressure drag" of kjto to get the 
total drag B due directly to the viscous forces. This 
strongly discredits the opinion sometimes advanced (as 
a result of unreliable tests at low Reynolds Number) 
that it is possible to get a negative pressure drag. It 
also makes the direct viscous drag of airships (and 
boats) quite sensitive to changes of form as well as to 
changes of surface area and ship speed. 

True viscous force. — To get an idea of the mag- 
nitude of force involved, O must now be evaluated, 
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taking into account the errors already introduced by 
the preliminary assumptions as to the exponents of 
the length and velocity. Using the formula of 
Diehl ("Engineering Aerodynamics," and N. A. C. A. 
Technical Note No. 102), 



0.01 



(7=0=0.0375 (6,350 X VJj)-° m ~ (y'JLf M 



(4) 



(for standard sea-level conditions); where L may be 
taken as the axial length of the hull, rather than the 
meridian length, to allow for the end taper; or, 

i=2a 

V, is the mean effective velocity, which may be 
evaluated as follows: 

Referring to equation (3), this may be written in the 
form, 



i2 = C7f SV. 2 



(30 
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Fineness ratio 
Fiotoe 12.— Variation of resistance with fineness ratio 

in which the surface area, 

S=~K, a b 

K, being a geometrical coefficient varying from 12.0 
for »=2 to 11.4 for n =8. From equations (3) and 
(3') we then have, 

K = C\a b K. F/=G|a s A* 

from which, 

T/2 T?A z n (n+2) V 0 * 
K * = K, (n+l) s 

Substituting in (4), 

(7=0.0077 



(n +2)° 



0.075 



(71+ l) 1 



and from (3), 

5 = 0.00011 A 3 -™ a 1 - 85 V 0 1Xi , 

71 

(for p=0.00238 and £, = 11.7) 

In terms of the volume (j* ^ g n i y' 

n (n+i) 0 - 16 

100A o:m vol 0 - 05 7l 0J7fi (7i+2)°- 07S V Q 0 - 16 



(5) 



(5') 



and, 



U ~ 22,000 (7i+ W 

For practical use this can be simplified to 

E = <7 s vol 0 - 817 V a 1M (6') 

where G B is a coefficient involving A and n. C<4 also 
depending on n) which is plotted on the attached 
chart (fig. 12) for various values of n. 

Discussion of results. — An interesting fact about the 
curve is that the minimuin drag is shown at a fineness 
ratio greater than what would cause any great amount 
of burbling over a properly shaped hull. Thus the 
fineness ratio for minimum drag is apparently deter- 
mined mainly by frictional forces alone. A supple- 
mentary indication is that practically minimum fric- 
tional drag, for a given volume, may be attained with 
any fineness ratio between about 3 and 6. But natu- 
rally the proper hull curve to prevent burbling becomes 
increasingly important for the lower fineness ratios; 
and with a ratio much under 4 a considerable amount 
of burbling must apparently be reckoned with at the 
best. In line with modern boundary layer theory, it 
appears further that this burbling at the lower fineness 
ratios is mdirectly due to the sharp increase of energy 
loss from viscous friction, particularly around the 
maximum section. 

The magnitude of the axial skin friction component 
around any one section, per unit axial length, is given 
by, 



^=2(7 2 ^LVcos 2 « 



(7) 



where O has the value given in equation (6). The 
approximation here involved in assuming O constant 

di? 

over a long hull makes the calculated value of a 

little too low at the bow and too high at the stern. 
Greater accuracy in respect to the actual force distri- 
bution will depend on further research into the detailed 
structure and growth of the boundary layer. 

Equations (6') and (7) are in a form which can be 
applied to any fair hull shape by estimating the equiva- 
lent ellipsoid, discussed in Part II. 

A similar method may be followed in dealing with 
other problems involving the boundary layer, includ- 
ing cases of mclined and curvilinear flight. 

It should be noted that there is still some question 
about the exact coefficient and exponent of the Rey- 
nolds Number. For smooth surfaces the coefficient 
as here used is more likely to be high than low. For 
the exponent, Von Karman uses — 0.2 instead of — 0.15 
(with a coefficient to correspond), but the bulk of the 
evidence seems to favor the exponent here used. (See 
reference 12.) 

Additional light on this point may be had from the 
values for various model tests in the variable density 



APPLICATION OF PRACTICAL HYDRODYNAMICS TO AIRSHIP DESIGN 



139 



wind tunnel as plotted in Figure 12. These are on the 
basis of Ob in equation (6'), all at a Reynolds Number 
of about 2X 10 7 . It will be noted that most of them 
fall below the theoretical curve; and most of them also 
show less variation with Reynolds Number than pre- 
dicted from theory. An offhand conclusion might be 
that the exponent and coefficient both need revision. 
However, data collected by B. M. Jones (reference 12), 
shows, quite clearly that the laminar type of flow is 
still often an appreciable element at the Reynolds 
Number of 10 7 , but probably would not be so at full 
scale values of 10 8 and over. In other words, these 
test results are probably not yet out of the transition 
stage, and hence do not quite comply with the first 
assumption on which this study was based. An experi- 
mental check of this supposition might be had from 
observing the change due to an artifiical increase of 
turbulence in the same tunnel. Lacking evidence to 
the contrary, it seems to be another case of the theoreti- 
cal result being on the whole more reliable and practical 
than the direct experiment; although in this case the 
extrapolation of several of the test results leads almost 
exactly to the theoretical value for full scale conditions. 

In the computation of the parasite drag of outside 
parts the proper coefficient and exponent will, of 
course, depend on the character of the part in ques- 
tion. If the true local velocity V (instead of V 0 ) is 
used, however, and the result further increased by the 
y 

energy factor ^ much of the additional drag at low 

fineness ratios, commonly attributed to interference, 
will be found accounted for. 

End conditions. — It remains finally to consider the 
instructive end conditions presented by the two ex- 
tremes of fineness ratio, 0 and <= . 

In the first case of n«=0, a circular disk normal to 

2tt 2 
the flow, it can be shown that K s a — and A= — > 

while L may be assumed equal to b (the radius). 
Substituting the corresponding value of G in equation 
(1) gives for the direct skin friction: 

i2o=(a finite coefficient) Xw 0 - 076 6 1 - 85 y 0 1 - 85 =0. This 
is not surprising in view of the fact that cos a 
is everywhere zero except at the extreme edge. But 
as energy is taken from the air, the drag can not be 
zero. On the energy basis, therefore, the viscous 
drag, here entirely in the form of "pressure drag," is, 

R=AR 0 = » 

a result involving flow dimensions obviously out of 
line with the original assumptions on which it is based, 
and hence inapplicable to a practical case. 

In the case of n= <», K b =t > , A -=1.0, and by a 
similar process, 

E=i2 0 = 0.0001l6a°- 8S To 1 - 85 
= 0.01 Sa (LV)- 0 - 16 



which leads back, as it should, to the basic skin fric- 
tion formula. 

CONCLUSIONS 

In conclusion of this part of the report, it is recom- 
mended for practical analysis that airship drag be 
divided into the following parts: 

1. Viscous hull drag, of which the proportion 

appears as skin friction proper, the balance as pres- 
sure drag. 

2. Burbling hull drag, creating an unknown propor- 
tion of pressure drag, usually positive, and a negative 
(relative to Item 1) skin friction drag. This item as a 
whole is probably always positive for practical hull 
shapes, but vanishingly small for the larger fineness 
ratios. 

3. Parasite drag of adjacent parts, considered as 
acted upon by the actual air flow locally in which 
they are placed. (See discussion above.) 

4. Interference effect, of the same parasite parts, on 
the hull drag (mostly modifying Item 2). 

Item 1 is obtained for full scale by formula (6') 
(the coefficient and exponents being subject to possible 
future refinement); or, for a model with unknown 
type of flow, by subtracting the axial resultant of the 
measured pressure distribution from the total meas- 
ured drag and multiplying the remainder by A. 

Item 3 is obtained by wind tunnel test or compu- 
tation corrected to speeds obtained by formula or 
estimate. 

The sum of Items 1 and 2 is obtained from a large 
or high-density tunnel, with maximum turbulence in 
the air flow. 

The sum of Items 1, 2, 3, and 4 is similarly obtained; 
and can also be had from actual flight tests. 

If C R is used as a coefficient of total drag, it has the 
following relation to the shape coefficient C s : 

, / _\ OJBS / „\ 0.16 0X17 

Drag = C s 2 vol 2 * = C>(j^) (±J vol. V 0 lM 

where po and ^ are the standard density and viscosity, 
respectively. 
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TABLE I 



INERTIA COEFFICIENTS OF ELLIPSOIDS 



a 




B— 

l+*» 


C= 
1+f 
«»+l 
n'-l 


B-A~ 
h-Ai 


AXB 


LO 


1.6 


1.6 


L0 


0 


2.25 


LB 


1.305 


1.621 


L244 


.316 


2.116 


2.0 


1.209 


1.702 


L400 


.493 


2.058 


2.5 


1.156 


L763 


L606 


.607 


2.038 


3.0 


1.122 


1.803 


L682 


.681 


2.023 


4.0 


L0S2 


1.860 


L689 


.778 


2.013 


6.0 


1.059 


L896 


L760 


.836 


2.007 


6.0 


L045 


1.918 


LS07 


.873 


2.004 


7.0 


L038 


L933 


1.839 


.897 


2.003 


8.0 


1.029 


L946 


L867 


.916 


2.001 


9.0 


L024 


L954 


L887 


.930 


2.001 


10.0 


L021 


1.960 


1.901 


.939 


2.001 


CO 


L0O0 


2.000 


2.000 


1.000 


2.000 



